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Using the self-consistent Hartree-Fock approximation for spinless elec- 
trons at zero temperature, we study tunneling of the interacting electron 
gas through a single d barrier in a finite one-dimensional (ID) wire con- 
nected to contacts. Our results exhibit features known from correlated 
many-body models. In particular, the conductance decays with the wire 
length as oc where the power a is universal. We also show that a 

similar result for a wire conductance can be extracted from the persistent 
current (I) through the 6 barrier in a ID ring, where it is known that 
I oc 

PACS numbers: 73.23.-b, 73.61. Ey 

Electron gas in a quantum wire is a realistic one-dimensional (ID) system. If 
a clean wire is biased by two macroscopic contacts with negligible electron back- 
scattering, the wire conductance is quantized as an integer multiple of e'^/h. This 
effect can be explained in the model of noninteracting electrons fL\. 

If a localized scatterer is introduced into the wire, quantization of conduc- 
tance breaks down due to the electron backscattering from the scatterer. For 
non-interacting electrons, the conductance is given by Landauer formula as the 
electron transmission probability at the Fermi level p]. However, the electron- 
electron (e-e) interaction alters the properties of the system qualitatively. In the 
Luttinger liquid model [2j, the conductance of an infinite wire containing a single 
scatterer varies with temperature as oc T~^" for T — > 0, where the power a depends 
only on the e-e interaction. For repulsive interaction a is positive and reflection 
at zero temperature is perfect, no matter how strong or weak the scatterer is. 

Matveev et al. [3] studied the Landauer conductance of the interacting ID 
electron gas through a 6 barrier in a wire with contacts. They replaced the many- 
body wave function by the Slater determinant of single-electron wave functions and 
analyzed the effect of the Hartree-Fock potential on the tunneling transmission. 
Assuming a weak e-e interaction of finite range, they derived the transmission 



using the renormalization group (RG). They confirmed the universal power law 
T~^", this approach is believed to go beyond the Hartree-Fock approximation. 

Here we consider the non-Luttinger liquid model of Matveev et al. [3j, but in- 
stead of the RG approach we apply the self-consistent Hartree-Fock solution. We 
evaluate the Landauer conductance. We find a good agreement with the theory of 
Matveev et al. In particular, we simulate asymptotic dependence of the conduc- 
tance on the wire length (L) for strong S barriers and we reproduce the universal 
power law oc We also show that essentially the same wire conductance can 

be extracted from the persistent current (J) in a ID ring, where / oc L~"~^. 

We consider a ID gas of interacting spinless electrons in a wire of length L. 
The wire is positioned along the x axis between x = —L/2 and x = L/2, both wire 
ends are connected to contacts. The single-electron wave functions ipk{x), where 
k is the electron wave vector, are described by the Hartree-Fock equation 



where j6{x) mimics the localized scatterer positioned in the center of the wire. 
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is the Hartree potential induced by the 6 barrier, 
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is the Fock nonlocal exchange term, and V{x — x') is the e-e interaction. 

The Landauer conductance is defined for the wire connected to large contacts 
via adiabatically tapered non-reflecting connectors [Tj. First assume a clean non- 
interacting wire, i.e., 7 = and V{x — x') = 0. As there is no backscattering at the 
wire ends, the solution of the equation ([1]) is the free wave e^^^ with eigenenergy 
h'^k'^/2m. The states e*'^^ and e"*'^^ with A; > describe the ballistic electrons 
originating from the contact at x = —L/2 and x = L/2, respectively. As they are 
mutually incoherent, e^^^ and e~*^^ are the only independent solutions. 

Second, let us keep 7 = but let us turn on the e-e interaction V{x — x'). 
Assuming that there is no backscattering at the wire ends, the solution of equation 
([T]) is still the free wave, il^^{x) = e^^^, but with the eigenenergy 



= h'e/2m + Ul{k), 
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where Up{k) = Uplipkix) = ipl{x)] is the Fock shift. Note that this solution is vahd 
if we imphcitly assume that the Fock interaction is present also in the contacts. 
Indeed, if the energy (jlj) holds inside the wire and we turn off the Fock shift to zero 
outside the wire, we obtain at each wire end the potential drop Up{k). This would 
cause backscattering at both wire ends and the solutions e'*^^ and e"**^^ would be 
no longer valid, in contrast with the ballistic conductance of clean wires [T]. 

Finally, consider the barrier ^S{x) and interaction V{x — x'). The 6 barrier 
induces Friedel oscillations of the Hartree-Fock potential. The electrons are thus 
scattered by the 6 barrier and by the oscillating potential relief. Since the scattering 
is elastic, the eigenenergy (jlj) remains unchanged and the wave function ipk{x) can 
be found by solving equation (JTj) as a tunneling problem with boundary conditions 

i^k{x = -L/2) = e'"^ + ne-''^, = L/2) = 4e*^^ (5) 

= -L/2) = t'^e-''^", = L/2) = e"*'^" + r^e''^", (6) 

where > 0, is the reflection amplitude, and tk is the transmission amplitude. 
Once we know the transmission, we know the Landauer conductance (e^/h) \tkpf- 

In reality the Friedel oscillations penetrate through the wire ends into the con- 
tacts, where they decay fast due to the enhanced dimensionality and decoherence. 
To mimic this decay within our ID model, we sharply turn off the oscillations 
to zero at both wire ends and keep Uh = and Up = Up{k) outside the wire. 
This constant ID potential emulates the non-refiecting connectors and justifies 
the above boundary conditions. Essentially the same ID model was a starting 
point of the RG study by Matveev et al. [3]. 

Unlike to Matveev et al., we solve the equation (JTj) by means of the self- 
consistent iterative procedure. To save the computational time and memory, we 
follow Ref. [1] and simplify the equation (j3j) as 

L/2 kp 

dx'V{x-x') / —Re[^l,{x')^k'{x)] (7) 
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by noticing that dk'ipy{x')ipk'{x) ~ 2t:6{x — x'). Unlike the exact form (jHj), 
the Fock potential (jTj) is local and independent on k. This saves time and allows 
us to simulate long wires. We present numerical results for the GaAs wire with 
electron density n = 5 x 10^ m~^, effective mass m = 0.067 mo, and e-e interaction 

V(x-x') = K)e-l"-"'l/''. (8) 

We adopt the finite-ranged interaction (jHj) because of comparison with the RG 
theory of Ref. [3] which also assumes the e-e interaction of finite range. Physical 
meaning of the finite range is screening. 
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Figure 1: Transmission probability at the Fermi level versus wire length for various 5 
barriers \tkp\ and various e-e interaction strengths a. The dashed curves show the data 
due to the RG formula ([9]); these data follow the asymptotic power law \tkp\ oc L~^° 
as the 5 barrier becomes strong. The circles connected by the full lines are our self- 
consistent Hartree-Fock data; they approach the same asymptotic power law. The powers 
a = 0.0277, 0.0561, and 0.0855 correspond to Vq = 11 meV, 22.3 meV, and 34 meV, 
respectively, with d fixed to 3 nm. If we adjust the same a by another choice of Vq and 
d, our results remain unchanged. 



According to Ref. j3j the bare 6 barrier, described by the transmission and 
reflection amplitudes tk and fk, is renormalized by the Friedel oscillations. The 
renormalized transmission probability at the Fermi level reads 



= ^fJ^f^'" - ^ {d/ir , (9) 

where d is the range of the e-e interaction V{x ~ x') and the right hand side of 
holds for small tkp and/or large L. For weak e-e interaction (a <^ 1) a reads 

VjO) - Vi2kj.) 
« = TTT ' (10) 

where V{q) is the Fourier transform of the e-e interaction V{x — x'). We evaluate 
a for our e-e interaction ([H]), for which V{q) = 2Vod/{l + q^dp). 

The bare amplitudes are tk = k/{k + i() and fk = + iQ, where C = 

^mjW'. Since kp and m are flxed, in the following we parametrize the bare 5 
barrier by its transmission coefficient \tkp\ ■ 

Figure [1] shows the transmission probability \tkp'^ versus the wire length L for 
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Figure 2: Self-consistent Hartree-Fock potential Uhf{x/L) along a ID wire, induced by 
the (not shown) 5 barrier at a; = 0. The 5 barrier is adjusted to have the transmission 
{ikp — 0-003 at the Fermi level (14 meV). The wire length, L, is varied as a parameter. 
The bottom panel shows just the positive maxima of the Friedel oscillations. 



various 6 barriers and various e-e interaction strengths. The RG formula © is 
presented by the dashed lines. For strong 6 barriers the dashed lines follow the 
asymptotic power law \tkpf oc in the log scale manifested by linear decay 

with slope —2a. Our Hartree-Fock curves (open symbols connected by full lines) 
show slightly higher transmission but clearly follow the same trend. In particular, 
for small enough \tkp\ and for not too small L all Hartree-Fock curves decay with 
the same slope as the RG curves, independently on the strength of the S barrier. 
It is instructive to show the Hartree-Fock potential Uh{x) + Uf{x) in the form 

UHF{y) = {d/Lr[UH{y) + UF{y)-U'^]/AL, y = x/L, AL = 7rhvF/L, (11) 

where we subtract the constant Fock shift Up in order to show exclusively the 
Hartree-Fock potential induced by the 6 barrier. 

In Fig. [2] we show the typical self-consistent Unpiy) in the wire with a strong 
scatterer at ?/ = 0. The potential exhibits Friedel oscillations with period \p/2L. 
The bare scatterer is thus "dressed" by an extra scatterer due to the Friedel oscil- 
lations. This is why we see the conductance to decay with L. It is more difficult to 
understand why we see just \tkp\ oc {d/Lf". As L increases, the Friedel oscilla- 
tions in Fig. [2] are too dense to be distinguishable, but we can observe asymptotic 
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Figure 3: (a) The IJhf data from the Fig El (6) Analogous data for the ID ring from 
our paper [5j. (c) Maxima of UuF for the wire (full line) and ring (dashed line), (d) 
The same as in the panel (c) but in log scale. 



decay of the oscillation amplitude. Notice that the "envelope" of the oscillation 
amplitude is the same for all L. Indeed, as shown in the bottom panel, the "en- 
velope" scales for all L to a single curve. Notice also (c.f. eq. (fTTl)). that Uupiy) 
involves the scaling factor [d/ L)^. This might be the reason why \tkp\^ oc [djVf''^ ^ 
but we have so far not found a clear interpretation. 

By tying the wire ends to each other one can create the ID ring. Magnetic flux 
applied through the opening of the ring gives rise to the equilibrium persistent 
current. In our paper [5] we have calculated the persistent current in the ID ring 
with a single 5 barrier. We have applied the same Hartree-Fock model as in this 
work, except that the Hartree-Fock equation with cyclic boundary condition [5] is 
the eigenvalue problem rather than the tunneling problem. 

In Fig. [3] we compare the Hartree-Fock potentials in the wire and ring. Obvi- 
ously, the amplitude of the Friedel oscillations in the ring saturates at the bound- 
aries. However, the amplitude of the Friedel oscillations in the wire decays with 
distance from the scatterer without any tendency to saturate (note in panel (d), 
that the decay of the full line in log scale is linear or even slightly faster than 
linear). Nevertheless, in both cases Uupiy) involves the scaling factor {d/L)'^ and 
the "envelope" scales to a single curve for all L. It is thus not surprising, that also 
the persistent current in the ring (/) comes out from the Hartree-Fock model [5j 
as a power law, namely LI oc {d/L)°'. This suggests that the conductance of the 
interacting wire might be obtainable from the persistent current in the interacting 
ring. In the non-interacting case / = (efi7'/2L)|tfc^| for magnetic flux 0.25/i/e. 
Applying this formula (intuitively) to the interacting electrons, we obtain the wire 
conductance from the formula jtfcj^P = 2LI /evp-, with / taken from the persistent 
current simulation [5] . In Fig. H] we show that the results (full symbols) reasonably 
agree with our direct calculation of the wire conductance. 
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Figure 4: The transmission probability for the ID wire (the open symbols represent 
the data from Fig. [1]) compared with the transmission probability extracted from the 
persistent current of the ID ring (full symbols). The left and right panels show the 
transmission probability for \tkp\ = 0.003 and \tkp\ = 0.139, respectively. 

In conclusion, using the self-consistent Hartree-Fock approximation at zero 
temperature, we have calculated the Landauer conductance of the weakly- interacting 
spinless electron gas in a ID wire with a single 6 barrier. We have found the uni- 
versal power law \tkpf oc [d/L)'^^ known from the Luttinger-liquid model [2] and 
RG models [31 E] ■ We have also found that essentially the same wire conductance 
can be extracted from the persistent current in a ID ring. 
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